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ABSTRACT 

We report on the first nulling interferometric observations with the Large 
Binocular Telescope Interferometer (LBTI), resolving the N’ band (9.81 - 
12.41 pm) emission around the nearby main-sequence star p Crv (F2V, 1-2 Gyr). 
The measured source null depth amounts to 4.40% ± 0.35% over a £eld-of-view 
of 140 mas in radius (~2.6 AU at the distance of rj Crv) and shows no significant 
variation over 35° of sky rotation. This relatively low null is unexpected given 
the total disk to star flux ratio measured by Spitzer/IRS (~23% across the N’ 
band), suggesting that a significant fraction of the dust lies within the central 
nulled response of the LBTI (79 mas or 1.4 AU). Modeling of the warm disk shows 
that it cannot resemble a scaled version of the Solar zodiacal cloud, unless it is 
almost perpendicular to the outer disk imaged by Herschel. It is more likely 
that the inner and outer disks are coplanar and the warm dust is located at a 
distance of 0.5-1.0AU, significantly closer than previously predicted by models 
of the IRS spectrum (~3AU). The predicted disk sizes can be reconciled if the 
warm disk is not centrosymmetric, or if the dust particles are dominated by very 
small grains. Both possibilities hint that a recent collision has produced much 
of the dust. Finally, we discuss the implications for the presence of dust at the 
distance where the insolation is the same as Earth’s (2.3 AU). 

Subject headings: circumstellar matter - infrared: stars- instrumentation; inter¬ 
ferometers - stars: individual [r] Crv) 


1. Introduction 

The possible presence of dust in the habitable zones of nearby main-sequence stars 
is considered a major threat for the direct imaging and characterization of Earth-like ex¬ 
trasolar planets (exo-Earths) with future dedicated space-based telescopes. Several inde¬ 
pendent studies have addressed this issue and concluded that visible to mid-infrared direct 
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detection of exo-Earths would be seriously hampered in the presence of dust disks 10 to 


(e.g., iReichman et ahN2006bl: iDefrere et al.N2010l: iRoberge et ahN2012[) . The prevalence of 
exozodiacal dust at such a level in the terrestrial planet region of nearby planetary systems 
is currently poorly constrained and must be determined to design these future space-based 
instruments. So far, only the bright end of the exoz odi luminosity function has been mea¬ 


sured on a statistically meaningful sample of stars (ILawler et al.l l2009l: iKennedv fc Wvatt 
2013h . Based on WISE observations and extrapolating over many orders of magnitude, 


Kennedy &: Wvatt fj2013 ') suggest that at least 10% of Gyr-old main-sequence stars may 
have sufficient exozodiacal dust to cause problems for future exo-Earth imaging missions. 
To determine the prevalence of exozodiacal dust at the faint end of the luminosity function, 
NASA has funded the Keck Interferometer Nuller (KIN) and the Large Binocular Telescope 
Interferometer (LBTI) to carry out surveys of nearby main sequence sta rs. Science observa¬ 


tions with the KIN started in 20 08 and the results were reported recently (IMillan-Gabet et ah 


20 111: iMennesson et ahllin pressi ). One of their analyses focused on a sample of 20 solar-type 
stars with no far infrared excess previously detected (i.e, no outer dust reservoir). Assum¬ 
ing a log-normal luminosity distribution, they derived the median level of exozodiacal dust 
around such stars to be below 60 times the solar value with high conhdence (95%). Yet, the 
state-of-the-art exozodi sensitivity achieved per object by the KIN is approximately one order 
of magnitude larger than that required to prepare future exo-Earth imaging instruments. 

The LBTI is the next step. A s urvey, called the Hunt for Obs ervable Signatures of Ter¬ 
restrial Planetary Systems lHOSTS. lHinal2013l:lDanchi et al.ll2014f) . will be carried out on 50 
to 60 carefully chosen nearby main-sequence stars over the next 4 years (Weinberger et ah in 
prep). In this paper, we report on the hrst mid-infrared nulling interferometric observations 
with the LBTI obtained on 2014 February 12 as part of LBTI’s commissioning. We observed 
the nearby main-sequence star rjCrv (F2V, 1.4 Mq, 18.3pc), previously known to harbor 
unusually high levels of cir cumstellar dust. The infrared excess was hrst det ected with IRAS 


fIStencel fc Ba.ckmanlll99ll) and con hrmed later by Syitzer/MIPS aX 70 /rm flBeichman et ah 


2006al) . Sy itzer/IRS at 5 - 35 u m f Ghen et al.l l2006l: Ihisse et al.ll2012l) . and Herschel at 70 
- 500/im flDuchene et al.ll2014l) . Interestingly, the SED exhibits two distinct peaks which 
can be explained by the presence of two well separated dust belts. The cold outer one was 
hrst resolved in the submillimeter and found to lie at a distance of ~150AU from the star 


fjWvatt et al.ll2005h . Mid-infrared interferometric observations with V LTI/MIDI and the 


KIN resolved the inner belt to be located within 3AU of the host star (ISmith et al.l 12008 


20091: iMillan-Gabet et al.ll20 


ll), a conclusion su pported by blackbody models that place the 


warm dust near lAU (e.g., IWvatt et al.ll2005l ). A detailed analysis of the dust composi¬ 
tion in the inner belt using the Spitzer/IRS spectrum suggests a location nearer to 3AU 
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due to the presence of graii is as small as 1 u m that are warmer than blackbodies for the 
same stellocentric distance flLisse et al.ll2012l) . Most importantly, the mid-infrared excess 
emission is remarkably larg e and intriguing considering th e old age of the system (1-2 Gyr, 


Ibukivama fc Arimotoll2002l:iMaHik et a].ll2nn3l: IVicanll2012l) and all the aforementioned stud¬ 


ies agree that it cannot be explained by mere transport of grain s from the outer b elt (at 
least not without a contrived set of planetary system param eters, Bonsor et all 2012 h This 
rather supports a rare and more violent recent event (e.g., iGaspar et ahll2013l) and makes 
Tj Grv an ideal object to study catastrophic events like the Late Heav y Bombardment that 
might have happened in the early solar system flGomes et ahl 120051) or a recent collision 
between planetesimals (e.g., as postulated for BD-1-20307. ISong et al.ll2005[) . 


2. Observations and data reduction 


2.1. Instrumental setup 

The Large Binocular Telescope consists of a two 8.4-m aperture optical telescopes on a 
single ALT-AZ mount installed on Mount Graham in southeastern Arizona (at an elevation 
of 3192 meters) and operated by an international collaboration among institutions in the 
United States, Italy, and Germany flHill et al.l I20l4 IVeillet et al.l 120141) . Both telescopes 
are equipped with deformable secondary mirrors which ar e driven with the LBT’s adaptive 


optic s system to correct atmospheric turbulence at 1 kHz flEsposito et al.ll2010l: iBailev et ah 
20141) . Each deformable mirror uses 672 actuators that routinely correct 400 modes and 


provid e Strehl ratios exceeding 80%, 95%, and 99% at 1.6/xm, 3.8 pm, and 10 pm, respec¬ 
tively flEsposito et al.ll2012l: ISkemer et al.ll2014ll . The LBTI is an interferometric instrument 
designed to coherently combine the beams from the two 8.4-m primar y mirrors of the L BT 
for high-angular resolution imaging at infrared wavelengths (1.5-13 //.m. lHinz et aLll2012h . It 
is developed and operated by the University of Arizona an d based on the h eritage of the 
Bracewell Infrared Nulling Gryostat on the MMT (BEING, Hinz et ah 2000). The overall 
LBTI system architecture and performance will be presented in full detail in a forthcoming 
publication (Hinz et al. in prep). In brief, the LBTI consists of a universal beam combiner 
(UBG) located at the bent center Gregorian focal station and a cryogenic Nulling Infrared 
Gamera (NIG). The UBG provides a combined focal plane for the two LBT apertures while 
the precise overlapping of the beams is done in the NIG cryostat . Nulling interf erometry, 
a technique proposed 36 years ago to image extra-solar planets (IBracewelll 119781) . is used 
to suppress the stellar light and improve the dynamic range of the observations. The basic 
principle is to combine the beams in phase opposition in order to strongly reduce the on-axis 
stellar light while transmitting the flux of off-axis sources located at angular spacings which 
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are odd multiples of 0.5A/i? (where B = 14.4 m is the distance between the telescope cen¬ 
ters and A is the wavelength of observation). Beam combination is done in the pupil plane 
on a 50/50 beamsplitter which can be translated to equalize the pathlengths between the 
two sides of the interferometer. One output of the interferometer is reflected on a short- 
pass dichroic and focu sed on the NOMIC camera (Nulling Optimized Mid-Infrared Camera, 


Hoffmann et al.ll2014j) . NOMIC uses a 1024x1024 Raytheon Aquarius detector split into 2 


columns of 8 contiguous channels. The optics provides a field of view (FOV) of 12 arcsecs 
with a plate-scale of 0.018 arcsecs. Tip/tilt and phase variations between the LBT apertures 
are measured using a fast-readout (IHz) K-band PICNIC detector (PHASECam) which re¬ 
ceives the near-infrared light from both outputs of the interferometer. Closed-loop correction 
uses a fast pathlength corrector installed in the UBC (see more details in Defrere et ah 2014). 


2.2. Observations 


Nulling interferometric observations of rj Crv were obtained in the N’ band (9.81 - 
12.41 /im) on 2014 February 12 as part of LBTFs commissioning. The basic observing block 
(OB) consisted of one thousand frames, each having an integration time of 85 ms, for a 
total acquisition time of 110 seconds including camera overheads. The observing sequence 
was composed of several successive OBs at null, i.e. with the beams from both telescopes 
coherently overlapped in phase opposition, and one OB of photometric measurements with 
the beams separated on the detector. In order to estimate and subtract the mid-IR back¬ 
ground, the OBs at null were acquired in two telescope nod positions separated by 2'.'3 on 
the detector. We acquired 3 different observations of rj Crv interleaved with 4 observations of 
reference stars to measure and calibrate the instrumental null floor (CAL1-SCI-CAL2-SCI- 
CAL1-SCI-CAL3 sequence; see target and calibrator information in Tabled]). To minimize 
systematic errors, calibrator targets were chosen close to the science target, both in terms 
of sky position and magnitude, using the SearchCal software flBonneau et al.ll201l[) . The 7 
observations occurred over a period of approximately 3 hours around the meridian transit 
(see u-v plane covered by the science observations in Fig. dj) with relatively stable weather 
conditions. The seeing was ~1'.'4 for the first observation and 1'.'0-172 for the remainder 
observations. The adaptive optics systems (AO) were locked with 300 modes (left side) and 
400 modes (right side) at a frequency of 990 Hz over the duration of the observation^. Fringe 
tracking was carried out at IkHz using a K-band ima ge of pupil fringes and an approach 
equivalent to group delay tracking (jPefrere et ahl 120141) . The phase setpoint was optimized 


Mhe different number of modes used on each telescope has a negligible impact on the null depth (i.e., a 
few 10“^). Furthermore, this effect is constant during the night and calibrates out completely. 
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Table 1: Basic properties of rj Civ and its calibrators. 


ID 

HD 

RA-J2000 
[d m s] 

DEC-J2000 
[d m s] 

Type 

mv 

rriK 

[Jy] 

^'lD ± Icr 
[mas] 

Refs. 

rj Crv 

109085 

12 32 04 

-16 11 46 

F2V 

4.30 

3.37 

1.76 

0.819 ± 0.119 

A13 

CAL 1 

108522 

12 28 02 

-14 27 03 

K4HI 

6.80 

3.46 

1.68 

1.204 ± 0.016 

M05 

CAL 2 

107418 

12 20 56 

-13 33 57 

KOHI 

5.15 

2.83 

2.96 

1.335 ± 0.092 

Bll 

CAL 3 

109272 

12 33 34 

-12 49 49 

G8HI 

5.59 

3.60 

1.42 

0.907 ± 0.063 

Bll 


References. Coordinates, spectral types and V/K magnitudes from SIMBAD; N’-band flux densities 
computed by SED fit (Weinberger et al. in prep); L imb-darkened diameters and 1-cr uncertainties 
from [A13] Absil et al. ( 201,31 ). [M05] Merand et al. (2005), and [Bll] Bonneau et al. ( 2011 ). 


at least once per observation to take ont any optical path variation between the K-band, 
where the phase is measnred and tracked, and the N’-band, where the nnll depth is measnred 
(making sure that the phase maintained at K-band provided the best possible nulls at N’ 
band). The NOMIC detector was used with the small well depth for better sensitivity and 
in subframe mode (512x512 pixels) to reduce the camera overhead. 


2.3. Data reduction and calibration 


Data reduction and calibration were performed using the nodrs pipeline developed by 
the LBTI and HOSTS teams for the survey (Defrere et al. in prep). It converts raw NOMIC 
frames to calibrated null measurements in hve main steps: frame selection, background 
subtraction, stellar flux computation, null computation, and null calibration. Frame selection 
is done by removing the hrst 20 frames of each OB that are affected by a transient behavior 
of the NOMIC detector. Background subtraction is achieved by nod pairs, subtracting from 
each frame the median of frames recorded in the other nod position. E ach row in a given 
channel is then corrected for low frequency noise (iHoffmann et al.l 120141) by subtracting its 
sigma-clipped median, excluding the region around the star position. Remaining bad pixels 
were identihed using the local standard deviation (5-0" threshold) and replaced with the mean 
of the neighbor pixels. The flux computation is done in each frame by aperture photometry 
using an aperture radius of 8 pixels (140 mas or 2.6 AU at the distance of r/Crv), equivalent 
to the half width at half maximum (HWHM) of the single-aperture PSF at 11.1 pm. The 
background flux in the photometric aperture is estimated simultaneously in each frame using 
all pixels covering the same columns as the photometric region and located between the 
second minimum of the single-aperture PSF (34 pixels or 600 mas) and the channel edge. 
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Raw null measurements are then obtained by dividing the individual flux measurements at 
null by the total flux estimated from the photometric OB of the same observation (accounting 
for the 50/50 beamsplitter). This procedure is performed exactly the same way for rjCiv 
and the three calibrators. Finally, in order to correct for the differential background error 
on the null due to the different brightness of the stars, we further add to the null a small 
fraction of background flux measured in a nearby empty region of the detector. This fraction 
is computed to match the background error on the null between the different stars, assuming 
fully uncorrelated noise in the photometric aperture and the nearby empty region of the 
detector (see Appendix |^. To compute the final null value per OB, we first remove the 
outliers using a 5-cr threshold and then perform the weighted-average of the lowest 5% null 
measurements (the weight being defined as the inverse square of the uncertainty in the 
aperture photometry). The corresponding null uncertainty is computed by bootstrapping 
through the entire dataset (keeping the lowest 5% each time) and taking the 16 and 84% 
levels from the cumulative distribution as representative 1-a uncertainties. A description of 
the additional systematic errors considered is given in Section |3l The raw null measurements 
per OB for t]Ctv and its calibrators are shown in the top panel of Figure [H 


The instrumental null floor is estimated using the calibrator null measurements corrected 
for the finite extension of the stars. This correction is done using a linear limb-darkened 
model for the geometric stellar null flAbsil et al.l 120061) : 


A^star — 




( 1 ) 


where B is the interferometric baseline, 0 ld is the limb-darkened angular diameter of the 
photosphere (see Table [1]), A is the effective wavelength, and ux is the linear limb-darkening 
coefficient. Given the relatively short interferometric baseline of the LBTI, the geometric 
stellar null will generally be small in the mid-infrared and the effect of limb-darkening neg¬ 
ligible. Assuming ux = 0, the typical geometric null for our targets is ~ 10“® with an error 
bar of ~ 10“®, which is small compared to our measurement errors. The next step in the 
calibration is to compute the instrumental null floor at the time of the science observations 
and subt ract it from the raw null measurements of the science target (null contributions are 
additive, IHanot et al.ll201l[) . This can be estimated in various ways, using only bracketting 
calibrator measurements, a weighted combination of the calibrators, or by polynomial inter¬ 
polation of all calibrator measurements. In the present case, because r] Crv and its calibrators 
have been chosen close in magnitude and position on the sky, we assume that the instrument 
behaved consistently for the duration of the observations and use a single/constant to fit the 
instrumental null floor, as shown by the solid line in the top panel of Figure [H The only 
noticeable difference between rj Crv (F2V) and its calibrators (G8III to K4III) is the spectral 










type but the slope of the stellar flux across the N’ band is the same so there is no differential 
chromatic bias on the null measurement. 

The statistical uncertainty on a calibrated null measurement is computed as the quadratic 
sum of its own statistical uncertainty and the statistical uncertainty on the null floor mea¬ 
sured locally at any given time by using the 1-cr variation on the constant £t of the instru¬ 
mental null floor. The systematic uncertainty, on the other hand, is computed globally and 
is composed of two terms. The first one accounts for the uncertainty on the stellar angular 
diameters, taking the correlation between calibrators into account. As described above, it 
is negligible in this case. The second one is computed as the weighted standard deviation 
of all calibrator measurements and accounts for all other systematic uncertainties, which 
are harder to estimate on a single-OB basis (see Section [3]). The statistical and systematic 
uncertainties are uncorrelated and added quadratically to give the total uncertainty on a 
calibrated null measurement. The part of this uncertainty related to the the null floor is rep¬ 
resented by the dashed lines in the top panel of Figure [TJ The calibrated null measurements 
are shown in Section [5l 


3. Data analysis and resnlts 

Figured] shows that the null measurements obtained on the science target (red diamonds) 
are clearly above the instrumental null floor (blue squares), suggestive of resolved emission 
around r] Crv. In order to quantify the measured excess emission, the hrst step of the 
data analysis is to convert the calibrated null measurements into a single value. Given 
that the calibrated null measurements show no signihcant variation as a function of the 
baseline rotation, this is done by computing the weighted average of all calibrated null 
measurements and gives 4.40%±0.10%±0.34%, where the two uncertainty terms correspond 
to the statistical and the systematic uncertainties respectively. The statistical uncertainty 
decreases with the number of data points and is computed as follows: 

_ 1 

C^stat j ■; (2) 

\J y ^/^i,stat 

where cTj^stat is the statistical part of the uncertainty on the ith calibrated data point. The 
systematic uncertainty on the other hand has been computed globally and we make here the 
conservative assumption that it is fully correlated between all data points. The systematic 
uncertainty on the final calibrated null is hence given by the systematic uncertainty on the 
null floor, which is the same for all calibrated data points. Several instrumental imperfections 
contribute to this uncertainty. First, there is a background measurement bias between the 
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Fig. 1.— Left, sampling of the Fourier {u,v) plane obtained for rjCiv on 2014 February 12. 
Each blue point and its centrosymmetric counterpart represent the orientati on of the LBTI 


basel ine for a given OB. The orientation of the outer disk midplane (i.e., 116?3. lDuchene et ah 


2014J ) is represented by the black dashed line. The right hgure shows the corresponding raw 
null measurements per OB as a function of UT time (top panel). The blue squares show 
the calibrator measurements while the red diamonds represent the rj Crv measurements. The 
estimated instrumental null floor is represented by the solid black line and the corresponding 
1-a uncertainty by the dotted lines. The bottom panel shows the corresponding background 
error estimate measured in a nearby empty region of the detector. 


photometric aperture and the nearby regions used for simultaneous background measurement 
and subtraction. The amplitude of this bias depends on the the photometric stability (i.e., 
precipitable water vapor, temperature, clouds) and the nodding frequency. It was estimated 
using various empty regions of the detector located around the photometric aperture. The 
result is shown for one representative empty region in the bottom panel of Figure [T] and 
accounts for 0.18% of the systematic uncertainty. Another main contributor to the systematic 
uncertainty comes from the mean phas e setpoint used to track the fringes. Using a more 
advanced data reduction technique fi.e.. lHanot et al.ll201ll) on similar data sets, we estimate 
that this error can produce a null uncertainty as large as 0.2% between different OBs. Finally, 


a variable intensity mismatch between the two beams can also impair the null floor stability. 
In the present case, it was measured in each observation using the photometric OBs and found 
to be stable at the 1.5% level, corresponding to a null error of ~0.1%. Adding quadratically 
these 3 terms gives a systematic uncertainty of 0.29%, similar to the value found after data 
reduction (i.e., 0.34%). 
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The final step of the data analysis is to compute the fraction of the calibrated source 
null depth that actually comes from the circumstellar environment by subtracting the null 
depth expected from the stellar photosphere alone. Using a limb-darkened diameter of 0.819 
± 0.119 mas (see Table [T]) and Equation [T|, the stellar contribution to the total null depth is 
0.0014% ± 0.0004% which is negligible compared to the calibrated null excess. Therefore, 
adding quadratically the statistical and systematic uncertainties, the final source null depth 
detected around rj Crv is 4.40% ± 0.35%. To make sure that no dust emission is lying outside 
the photometric aperture, we reduced the data using larger apertures and following exactly 
the same procedure. The results are illustrated in Figure [2] for four different aperture radii 
which are multiple of the half width at half maximum of the single-aperture PSF at 11.1 /rm. 
The lack of significant increase in the null measurements for larger aperture radii suggests 
that the angular size of the inner disk is smaller than the size of the single-aperture PSF 


Smith et ah 

2008 

) and the KIN ( 

Mennesson et ah 


m pres: 


3 ). 


4. Modeling and interpretation 


A first step toward interpreting the calibrated null depth observed for r] Crv can be 
made simply by comparing it with the ratio of disk to phot o spheric flux observed with 
Spitzer/IRS. Using the spectrum presented in Lebreton et ah f in prepi ) and assuming an 
absolute calibration error of 3%, we derive a value of 22.7% ± 5.4% over the N’ band. 
This is a very conservative uncerta inty estimate which is consistent with previous studies 
fjChen et al.ll2006l : iLisse et al.ll2012h . For any centrosymmetric disk, approximately half the 
disk flux is transmitted through the LBTI fringe pattern unless the extent of the N’-band 
emission is comparable to, or smaller than the angular resolution of the interferometer (i.e., 
A/2i?=79mas or 1.4 AU). The observed null depth of 4.40% ± 0.35%, i.e. substantially 
lower than the disk to star flux ratio, suggests therefore that the disk is either compact 
or edge-on with a position angle perpendicular to the LBTI baseline. An edge-on disk 
perpendicular to the baseline would have a position angl e of around 0-20° w hich is specifically 
disfavored by VLTI/MIDI observations (see Fig. 11 in ISmith et al.ll2009h . In addition, the 
inner disk would be nearly perpendicular to the outer disk plane whose o rientation has been 
consistently measured by independent studies (see dashed line in Figure fll I Wvatt et al.ll2005 


Duchene et al.ll2014J: iLebreton et ah 


in predL Given that the w’arm dust may be scattered 


in from the cool outer belt (e.g., IWvatt et al.l 120071: iLisse et al.ll2012l) . such a configuration 


seems unlikely and would be difficult to explain dynamically. Besides, among the planetary 
systems with measured inner and outer disk orientations (e.g.. Solar system, /3Pic, AUMic), 
none shows perpendicular inner and outer disks. A more likely scenario is for the inner and 
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Photometric aperture radius iAU] 

0 2 4 h 8 10 12 



Fig. 2.— Calibrated null excess as a function of photometric aperture radius. The left 
vertical dashed line indicates the interferometric inner working angle (i.e., A/4i?=39mas or 
0.7 AU at the distance of rjCiv) while the right vertical dotted line represents the default 
aperture radius used in this study (i.e., the half width at half maximum of the single-aperture 
PSF at 11.1 fim). The lack of signihcant increase in the null measurements for larger aperture 
radii suggests that the regions emitting the majority of the N’-band flux is more compact 
than the size of the single-aperture PSF. 


outer disk components to be coplanar, and to concentrate more than half the disk flux within 
1.4AU. 


In order to interpret our observations within the context of a physical model, we con¬ 
struct a simple vertic ally thin model for the inner disk using the analytical tool developed 
by the HOSTS team flKennedv et ahl [submitted to ApJI) . As expected from the qualitative 
description above, scaled models of the Solar zodiacal cloud that match the observed total 
and transmitted disk fluxes have position angles near-perpendicular to the outer disk. That 
is, this model is too radially extended so, unless the inner disk is strongly misaligned with 
the outer disk, a scaled version of our Solar zodiacal cloud is not a good match to the warm 
disk around rj Crv. We therefore explore the alternative possibility that the inner and outer 
disks are coplanar. With the disk orientation fixed, it is necessary to vary the location and 
width of the warm dust disk to reproduce both the total disk flux and the calibrated null 
measured by LBTI. That is, the only free parameters are the disk radius and width (r and 
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Fig. 3.— Constraints on the location and width of the resolved excess emission detected 
aronnd rj Crv. Grey contonrs and labels show the predicted LBTI nnll assnming a disk 
to star flnx ratio of ~23% and the same orientation as the onter belt. Bine lines show 
the observed LBTI nnll and corresponding nncertainties. The left red line shows the inner 
disk radins lower limit of 0.16 AU as derived with VLTI/MIDI while the right dark red line 
indicates the onter limit set by onr observations. Considering all constraints, the disk is 
restricted to he at a distance of 0.5-1.0 AU, with a fnll width of np to 1.7 AU. 


Sr), and the power-law for the dnst (i.e., snrface density or optical depth) between these 
radii (a). The resnlts of this process, for the radins and width, are shown in FignreO The 
greyscale and grey contonrs show the nnll depths of models with a range of radii r and 
fractional widths Sr/r, where in each case the total disk to star flnx ratio is 23% in the N’ 
band (nsing a = 0.34). The nnll depths vary from near zero for small disk radii, to abont 
11% for narrow rings with r 1.5 AU. The bine lines show where the predicted nnll depths 
agree with the LBTI observations of rj Crv with the range allowed by 2 and 3-(T nncertainties. 
These nncertainties do not acconnt for the 5% nncertainty in the Spitzer/IRS excess, which 
for example moves the bine lines systematically approximately 0.1 AU ontward for a 5% 
lower excess. Similarly, if the IRS excess is lar ger, which is possible if the excess is present 
at shorter wavelengths (e.g., iLisse et aLN2012h . then the disk lies even closer to the star. 
With these assnmptions, the inner disk component is constrained to he at approximately 
0.5-1.0AU with varying width. 













13 



Fig. 4.— Surface density as a function of distance from the star for five disk models having 
different power-law indices a. Each line represents the widest disk allowed in a Figure [3] 
computed for the corresponding a. For a steeply decreasing surface density profile (i.e., 
large a), the disk can extend beyond 1AU but the disk emission is then more concentrated. 


As noted in Section [3l the lack of increase in the null measurements for larger aperture 
radii suggests that the d isk is smalle r than 2.6 AU in radius. From MIDI, VISIR, and 
MICHELLE observations. ISmith et al.l (1200911 constrain the inner disk to lie outside 0.16 AU. 
These constraints are shown by the red curves in Figure |3l and serve as useful limits if the 
disk is wide, in which case the LBTI does not constrain the disk extent strongly. The 
constraint is in fact better than appears from Figure [3] because even the widest models are 
concentrated; half of the total disk flux originates from within 0.8 AU. Wide disks are only 
allowed because enough of their emission is blocked by the LBTI transmission pattern so 
the qualitative description given above holds. The zoo of possible models can be made even 
larger by adjusting the surface density power law. The default value of a=0.34 used above 
comes from the solar Zodiacal cloud model flKelsall et al.l 1199811 which is not necessarily a 
good match to the warm disk around rj Crv. It could be very different depending on the 
physical origin of the disk, and depends on dynamical, collisional, and radiation processes. 
Given that the warm dust around p Crv may be scattered in from the cool outer belt (e.g., 
Wvatt et ah 2007 : Lisse et ^12012 1. predictions of the expected a would be highly uncertain. 


However, its value has a very weak impact on the warm dust location as shown in Figure 01 
For models where the surface density increases with radius {a < 0), the disk emission is 
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concentrated farther from the star, so the disk mnst be smaller in order to £t our data. 
If the disk surface density decreases with radius (a > 0), our ability to constrain the disk 
extent is limited. However, for these models most of the emission comes from the inner disk 
edge, so the inability to constrain the outer edge location is expected, and applies similarly 
to other observations. Therefore, regardless of the disk profile, the conclusion from any 
centrosymmetric model is that the bulk of the disk emission lies closer to the star than 
previously thought, most likely in the range 0.5-1 AU. 


5. Discussion 


The warm dust location derived in the previous section is signific antly closer to th e 
star than inferred from Spitzer/IRS spectrum models (3AU or 160 mas. iLisse et al.ll2012h . 
There is however significant uncertainty in this inference because it relies on complex grain 
models with degenerate parameters. While the LBTI observations clearly suggest that these 
models need to be revisited, an alternative scenario to reconcile the difference in the inferred 
disk sizes is to relax the assumption of a centrosymmetric disk. If the inner disk has an 
over-density that was hidden behind a transmission minimum for the range of observed 
hour angles, then the disk can be larger and may still satisfy the observed LBTI null and 
Spitzer/IRS excess (i.e., a clump can be at larger physical separation, but remain close to 
the star in sky-projected distance). Future observations can therefore aim to expand the 
range of hour angles and also look for significant variations in the calibrated null depth at 
the same hour angles. Such observations will help establish whether an asymmetry exists, 
and if so, whether that asymmetry is orbiting the star or fixed in space. The latter possibility 
would favor recent collision scenarios which suggest that a long- lasting asymmetry is created 
at the location where the impact occurred fjJackson et al.ll2014ll . 


New modeling based on Spitzer/IRS models and KIN observation s also su ggests a com¬ 
pact warm disk emitting mostly from wi thin lAU flLebreton et ahl lin prenh. Using the 
GraTer radiative transfer modeling code flAugereau et ah 1999 ). Lebreton et ah f in prep ) 
propose two inner disk models that reconcile the KIN observations with the shape of the 
mid-infrared Spitzer/IRS spectrum and various photometric constraints from IRAS, AKARI, 
WISE, and Syitzer/ MIPS. The first m odel (model 1 hereafter) consists of small amorphous 
forsterite fMgoSiO^. Hager et al.ll2003l) dust grains with a differential size distribution defined 
by n(a) oc a~'^, with k, = 4.5 and Omin = 1.2 pm. The density peak location of the grains is 
0.16-0.25 AU with a radial decrease in surface density oc Note that forsterite grains are 
also the dominant species in the Spitzer/IRS best-fit model of iLisse et ahl (120121). Th e sec¬ 
ond model (model 2 hereafter) consists of astronomical silicates flDraine fc Leelll984j) . with 









































15 


14 



0 10 20 30 

Parallactic angle [deg] 


Fig. 5.— Calibrated null measurements as a function of parallactic angle for our LBTI ob¬ 
servations of r] Crv (Fe bruary 12, 2014') . The nu ll expected from the two best-fit exozodiacal 
disk models derived by iLebreton et al.l flin prepi ) to reconcile the KIN obseryations with the 
shape of the Spitzer/IRS spectrum are shown by the grey-shaded areas (see main text for 
model description). The width of the grey-shaded areas is set by the uncertainty on the 
density peak location of the model. 


a minimum size of Omin = 0.6 /rm and a size distribution slope oi k = 3.5. The density peak 
location of the grains is 0.73-0.85 AU with an outer density slope of —1. The expected null 
excess for both models is shown in Figure [5] as a function of parallactic angle. As indicated 
by the grey-shaded areas, model 1 is a better fit to the LBTI observations, particularly for 
its larger value of the density peak position (0.25 AU). To understand how the LBTI ob¬ 
seryations break the model degeneracy, the two disk models are shown on top of the LBTI 
and KIN transmission maps in Figure [6l While both models logically produce a similar flux 
at the output of the KIN, the expected flux transmitted to the null output of the LBTI is 
signihcantly lower for model 1 due to the shorter nulling baseline and the small extent of 
the disk model. Approximately 10-19% of the disk flux is transmitted to the null output for 
model 1, depending on the density peak position (i.e., 0.16-0.25 AU), while this value reaches 
48-51% for model 2 (0.73-0.85 AU). This explains the difference in the predicted null excess 
shown in Figure O While the parameters of model 1 could be tuned to better fit our ob¬ 
servations, this analysis clearly shows that it is possible to reconcile the Spitzer/IRS model 
with spatial constraints from the KIN and the LBTI. The very steep size distribution and 
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resulting large population of very small grains of such a model depart signific antly from the 
behavior expected for a collisional cascade in equilibrium flGaspar et al.ll20121) . It is possible 
that this is evidence for t he dust being the product of a transient, high rate of collisions 
(e.g., Olofsson et ah 2012 1. 


An important question related to the overall goal of LBTI is what the dust level in the 
habitable zone around rj Crv is. From the above discussion, it is clear that a scaled version 
of the Solar zodiacal dust cloud is not a g ood match to the o bservations and quoting a 
“zodi” level could therefore be confusing (see iRoberge et al.ll2012l. for a discussion about this 
unit). To quantify the amount of warm dust, we can instead compare the surface density 
of our models with that of the Solar zodiacal cloud at the same equivalent position (i.e., 
scaled as -v/TT). Several possible scenarios have been discussed, and each will make different 
predictions. In the simplest scenario, where the inner disk is coplanar with the outer disk, 
the analytical model presented in Figure[3]has a surface density at 0.7 AU approximately 10"^ 
times greater than that at the equivalent location in the Solar zodiacal cloud (i.e., 0.3 AU). 
At 2.3 AU, the distance where the insolation is the same as Earth’s at ^/Ll, it has formally 
no dust since models within the 2-a contours in Figure [3] all have outer edges within 2AU. 
The radiative transfer disk model (model 1) on the other hand is not truncated at its outer 
edge (surface density oc r“^'®) and has a surface density at 2.3 AU approximately 12000 times 
larger than that at lAU in the Solar zodiacal cloud. Finally, in the non-centrosymmetric 
scenarios, the possibility of even more dust at larger radii means that the habitable zone 
dust levels could be even larger, but localized to certa in azimuthal ranges which would be a 


source of confusion for future imaging missions (e.g., IDefrere et al.ll2012l) . Even if the lack 


of variation of the null measurements as a function of the observed baseline position angle 
suggests that the disk is most likely centrosymmetric, this scenario cannot be ruled out with 
the present data and more observations expanding the range of sky rotation are required. 


6. Summary and conclusions 

This paper presents the first science results from the LBTI nuller, a mid-infrared in¬ 
terferometric instrument combining the two 8.4-m apertures of the LBT. As part of LBTI’s 
commissioning observations, we observed the nearby main-sequence star rj Crv, previously 
known to harbor high levels of circumstellar dust. The data are consistent with a source 
null depth of 4.40% ± 0.35% across the N’ band (9.81 - 12.41 pm) and over a field-of-view 
of 140 mas in radius (~2.6AU at the distance of rj Crv). The measured null shows no sig¬ 
nificant variation over 35° of sky rotation and is relatively low compared to the total disk 
to star flux ratio of 23% observed with Spitzer/IRS, suggesting that a significant fraction 
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Fig. 6.— The two best-fit exozodiacal disk models derived by iLebreton et al.l fjin prepi) to 
reconcile the KIN observations with the shape of the Spitzer/IRS spectrum shown on top 
of the sky transmission map of the LBTI (top row) and the KIN (bottom row). Model 1 is 
on the left and model 2 on the right (see main text for model descriptions). The contours 
mark regions corresponding to given ratios of the maximum disk flux density (see contour 
label), assum ing that the inner wa rm disk has the same orientation as the outer disk imaged 
by Herschel I Dnchene et ah 20141) . While the predicted flux transmitted through the KIN is 
logically similar for both models (see value in the upper left legend), it is signihcantly different 
in the case of the LBTI which is well matched for few AU-scale dust emission around near 
main-sequence stars and allows to break the degeneracy between the models (see Figure [5]for 
the corresponding expected null measurements). The LBTI sky transmission map is shown 
for the first OB of our observations (Julian date: 2456700.9204, parallactic angle: -4?9, hour 
angle: -0hl8). The KIN sky transmission map is shown for the same hour angle. North is 
up. East is to the left. 
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of the dust lies within the central nulled response of the LBTI (79 mas or 1.4 AU along the 
baseline orientation). Assuming that the inner and outer disk components are coplanar, we 
show that the bulk of the warm dust emission must lie at a distance of 0.5-1.0AU from the 
star in order to £t our data. This is consistent with spatial constraints from the KIN that 
also point toward a compact warm disk emitting mostly from within 1 AU but significantly 
closer than the distance predicted by the Spitzer/IRS spectrum models (3AU). This dis¬ 
crepancy illustrates how spatially resolved observations provide crucial information to break 
the degeneracy inherent to SED-based modeling of proto-planetary and debris disks showing 
evidence of warm emission . This is discussed in de tail for r] Crv in a companion paper based 
on the KIN observations f Lebreton et ah in prepi) . In this paper, we show how the LBTI 


observations break the degeneracy between models that have the same SEDs and KIN nulls 
and discuss an alternative scenario which is the possible presence of an over-density in the 
inner disk. Both scenarios support the prevailing interpretation for the origin of the warm 
dust which is a recent collision in the inner planetary system. To conclude, it is clear that 
Tj Crv is not a good target for a future exo-Earth imaging mission since the surface density 
of the dust in or near the habitable zone can be as high as four orders of magnitude larger 
than that in the Solar zodiacal cloud. 
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A. Correcting the background bias between stars of different magnitudes 


The best 5% approach used in this paper to convert a series of individual null measure¬ 
ments to a single value creates a bias between stars of different magnitudes. The origin of 
this bias is the background noise which is presumably constant in absolute ffux and divided 
by the ffux of the star to compute the null. Therefore, the width of the relative background 
ffux distribution decreases with the brightness of the star and the mean of the best 5% mea¬ 
surements increases. The contribution of the background noise to the measured null then 
creates a bias that depends on the brightness of the star. In order to compensate for this 
effect, we can add a small fraction of background ffux to the measured ffux at null to match 
the relative background noise on stars of different magnitudes. Explicitly, we are looking for 
the fraction x verifying the following equation; 


variance 



= variance 



J- xB{t) 


(Al) 


where Fi and F 2 are the stellar ffuxes measured in the photometric OBs (assuming Fi > F 2 ), 
Bi{t) and B 2 {t) are the instantaneous background ffuxes in the photometric aperture in the 
null OBs, and B{t) is the instantaneous background ffux measured in a nearby empty region 
of the detector. Assuming that Bi(t) and B(t) are uncorrelated, their variance adds up and 
the above equation becomes: 


1 



X = 


F1F2 


(A2) 


This preprint was prepared with the AAS lATRX macros v5.2. 
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where we have assumed that (Jb^= cjb 2 = <^b- 



